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Abstract 


The assumptions of a region of slip near a moving contact line (to remove the force* 
singularity) and a constant contact le are used to obtain the equation for the shape of 

a thin drop of liquid resting on a plane. Three asymptotic expansions are matched 

together to obtain an expression foi cn«. rate at which the drop spreads. Some cases of 
sliding motion are also examined. Alv: ugh the technique is presented here for thin drops 
only, it can also be applied to drops of arbitrary size. 

Introduction 


The chief difficulty that has to be overcome in an att^m^t to describe the motion of a 
drop in contact with a rigid surface is the determination of the correct boundary condi* 
tions to be app.i ied at the contact line where the surface of the drop meets the rigid sur- 
face. If the drop is spreading over the surface, or if it is sliding along it, or both, 
the contact line is moving and it is well known that in these circumstances application of 
the no-slip boundary condition leads to a solution containing an unacceptable singularity. 
For certain purposes this singularity can be ignored but if, for example, we wish to deter- 
mine the rate at which the drop spreads or the speed at which it slides, a dynamical bal- 
ance of the forces acting must be achieved, which proves to be impossible in the presence 
of the force-singularity. To circumvent this difficulty, the most widely used device is to 
replace the no-slip boundary condition in the vicinity of the contact line by one allowing 
a certain amount of slip there. The argument in favour of this proposal is that large 
stresses occur near the contact line, associated with the rapid change in direction of the 
fluid motion there, and the molecular forces which are usually sufficient to prevent aay 
slip between fluid and solid may be unable to control these large stresses. Of course, 
what is really required is a good molecular theory for the junction between drop, solid and 
surrounding air. Failing such a theory, some model boundary condition may be proposed in 
the hope that its exact form is of little consequence and that it will at least enable 
finite answers to be produced in answer to questions involving the force balances associ- 
ated with the motion and that these answers can then be tested experimentally. The sim- 
plest model boundary condition is one which replaces the usual no-slip condition by one 
which allows a small amount of slip, proportional to the local velocity gradient. If u is 
the velocity parallel to the plane surface and z is measured normally away from it, the 
proposed boundary condition is 


u -13 u/3z =0 (1) 

for a fluid in contact with a solid boundary at rest. The slip coefficient 1 is a measure 
of the length over which slip is significant. Although its value is unknown, if the slip 
is produced by molecular effects we might expect a value of the order of 10”^ m. Such 
slip is only expected to occur close to the contact line and the boundary condition should 
revert to its usual form elsewhere. However, the very small size suggested for X indicates 
that the slip condition can be used everywhere with negligibly small error. As we shall 
see, allowing for slip near the contact line gives speeds proportional to l/|inX| which 
is of much greater significance than any erroneous inclusion of terms proportional to X . 

When we have decided to use a slip boundary condition for problems involving moving con- 
tact lines, the stresses near the contact line are still large and we may consequently 
expect that there will be a significant distortion of the shape of the free surface near 
the contact line. The large stresses will be balanced by capillary effects, which are pro- 
portional to the curvature of the surface. In order to provide sufficient boundary condi- 
tions for the shape of the free surface to be calculated, it is necessary to specify the 
angle at which the free surface meets the plane* Observations of contact angles seem to 
show that the angle increases with the speed of advance of the contact line, and decreases 
wi.^n the contact line is retreating. For stationary drops, the contact angle is not 
uniquely defined but can take any value between certain limits. The evidence for the 
dynaonic behaviour of the contact angle is based on observations which do not take account 
of any rapid change of slope in the very small distance from the edge of the drop where, 
as we have seen, large stresses are present. It may well be that the contact angle 
measured at the edge itself does not change with speed and the dynamic behaviour refers 
only to an apparent contact angle, only relevant at some distance from the edge. Evidence 
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in favour of this contention has been provided by Lowndes^ who has produced nuxoerical cal- 
culations for the motion of a meniscus along a tube using the slip hypothesis and keeping 
the actual contact angle fixed. He was able to show that apparent contact angles could 
be found from his calculated meniscus shapes which were in good agreement with observed 
values and that marked changes in slope occurred in the iiwnediate vicinity of the edge. 

The calculations of Lowndes ^ show that a self-consistent rational framework for dealing 
with moving contact lines is provided by the assumptions of slip at the edge and a fixed 
contact angle. This fr 2 unework has been used to discuss the spreading and sliding of a drop 
which is thin enough for lubrication theory to be used^. These problems are time-dependent 
and an estimate of the speeds involved depends crucially on the slip hypothesis. Both the 
meniscus problem and the drop problems needed considerable numerical calculation, although 
the use of lubrication theory reduced the amount required to a large extent. The equations 
for the spreading and sliding drop problems involve a small parameter, the slip coefficient, 
which suggests that matched asymptotic expansions could yield the desired answers without 
recourse to extensive numerical calculation. An attempt at such a procedure was made 
before the numerical results were obtained, but was not successful. The present paper, 
however, shows that matched asymptotic ea^ansions can be used to obtain the rate of spread 
of a thin drop and, with certain restrictions, the rate at which a drop slides down a plane. 
The key to success was the realisation that, as well as expected inner and outer expansions, 
an intermediate region was required across which the inner and outer regions could be 
matched. The method employed has some similarities to the one set out by Lacey^, except 
that he used a multiple scale approach and did not carry the solution sufficiently far. 

Although the methods to be described here can be applied to a variety of problems, for 
simplicity they are explained with reference to the problems examined before^. The drop is 
assumed to be thin and two-dimensional and it is also assumed that it is small enough for 
the Bond number to be small, that is, for gravity to be less significant than capillarity. 

In the first problem, the drop is placed on a horizontal plane and allowed to spread until 
its equilibrium position is reached. The quantity to be determined is the rate of spread 
of the drop as a function of its width. In the second problem, the drop is placed on an 
inclined plane, when both spreading and sliding may occur. The final width of the drop and 
the speed at which it slides are the quantities to be determined in this case. 

Formulation 


The application of lubrication theory to the problems to be solved is straightforward 
and has been described in detail before^. The simplified forms of the Stokes equations 
enable the velocity components to be found in terms of the pressure and the application of 
the normal stress condition and the kinematic boundary condition at the free surface yields 
an equation for the height h(x,t) of the drop,. where x is measured parallel to the plane 
surface and t is the time. The equation for h is 

O at 3x L O ;j 

where y is the fluid viscosity, g gravity, a surface tension, X the slip coefficient 
and 0 the inclination of the plane to the horizontal. Lubrication theory is only valid 
if the slope of the drop surface is everywhere small, so the contact angles must be small. 

A non-dimensional form of (2) is 

[ h2(h. X) j0-K)] - 0 , (3) 

;*;here 

h =* a^o h, x ■ ax, t ** 3wa^t/oa^, X • a a^X/3, K *■ a?pg sin 6/oo^ , (4) 

oa o oa oa o a 

and a is a length scale associated with the size of the drop and a is a typical value 
of the^ (small) contact angle. 


The extent of the plane covered by the drop can be fixed by two more unknowns, a. (t) 
and a^Ct), so that a 2 4 x < a^ and 

h(a^(t) ,t) - ,t) “ 0. (5) 

The volume of fluid in the drop remains constant throughout the motion and the length scale 
can be chosen so that 


a j 

h\x,t) dx ■ ^ 


(6) 
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The other conditions on h relate to the slope of the surface of the drop at its edges. 

The proposal indicated in the Introduction shows thatr when thet contact line is moving in a 
direction from th interior of the drop to the exterior (an advancing edge) , the contact 
angle is fixed at its static value, which we can take to be the scaling factor a . For 
a retracting edge, when the motion is in the opposite direction, the contact angle is equal 
to the minimum static angle . At a stationary edge the contact angle can lie anywhere 
between and . The boundary conditions to be applied at the edges are, therefore, 


at X = a, : - 3h/3x = 1 

B ^ - 3h/3x ( 1 

“ 3h/3x = B 

at X * a«: 3h/3x « B 

B ^ 3h/3x ( 1 

3h/3x = I 

where B = and 0 ^ B ^ 1 . 

1 . a 


if 

da. /dt 

> 

0, 


if 

dat/dt 

s 

0, 

(7) 

if 

daj/dt 

< 

0, 


if 

da-/dt 

> 

0, 


if 

daVdt 

s: 

0, 

(8) 

if 

da‘/dt 

< 

0, 



The final information needed to specify the problem completely is the initial state, 
that is, the values of a.(0), a2(0) and h(x,0) . From arbitrary initial states, one 
expects a fairly rapid transient'^phase during which the shape of the drop changes without 


significant spreading taking place. This is because the rate of spread is controlled by 
conditions at the contact line where large stresses resist the motion. No such restriction 
is placed on the distortion of the drop surface when the edges are fixed. Since the 
initial phase is of little significance to the general spreading problem it will be ignored 
here. The height of the drop is then a function of time only through its dependence on the 
positions of the edges. Thus the solutions to be obtained are of similarity type and the 
initial shape of the drop will be assumed to agree with these solutions for the given 
initial positions of the edges. 


Spread i ng 


The simplest problem of the type being considered is when the drop spreads on a 
horizontal surface. If we suppose that the drop starts from a position wher^ the apparent 
contact angles are greater than the static value, the edges of the drop will move outwards. 
If the drop is initially symmetric about its raid-position, it will remain so and we may 
write a. (t) = -a 2 (t) = a(t) and consider the interval 0 jp x < a only. The problem then 
takes the form 


• Is * w * '' 


(9) 

h = 0, -3h/3x *-- 1 

at x = a , 

(10) 

3h/3x = 0 

at X = 0 , 

(11) 

a 



h dx = 1/3, 


(12) 


o 


where a = da/dt and h is a function of x and a , in line with the intention to ignore 
any initial transient phase. As already explained, the rate of spread is expected to be 
small so that we caffi expand in powers of some small parameter c(X) , where X << e << 1 
but the dependence of c on X is yet to be determined. Hence we can write 


a = €Uj + ' 


(13) 


and, in the outer region where x is 0(1), h *» h^ + ch^ + ... 
equations satisfied by h^ and hj are 


h[-o 


- 

3x3 . 

3x3 


0 , 


3a 


From equation (9) the 


(14) 

(15) 


and these have to be solved subject to the conditions (10), (11) and (12), except that the 
condition on the slope at the edge of the drop is not to be applied in this outer region. 
The solutions are 


$ 


\ 
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^1 “ ^1*** tn(a + x) + (a - x) 4 n(a - x) - 2a in 2a + 


so that, near the edge of the drop where a * x is small. 


h 'v, - ~ 2 ^ ■♦■ c u^a**(a + x) fin ^5^ + 2 ] . 


The inner region is close to the edge of the drop and is where slip is important. In 
this region we write 

X * a - AX, h(x,a) - AH(X,a), H = Hq + cH^ + ... { 

The equation satisfied by H witn a term 0(A) omitted is 


^ 37 [h^‘h ^ If] = ° 


and the boundary conditions are 

H « 0, 3H/3X « 1 at X * 0 , (21) 

and 3H/3X is not more than logarithmically infinite as X-^« . The solution is 

H = X + ijeu^ (X + 1)2 in(X + 1 ) - X^ tn X - X) j + . . . , (22) 

SO that, toy’ X-*-« , 

H -x. X + £U^X in X . (23) 

Although (l8) and (23) both contain logarithmic terms, it is not possible to match these 
two expansions together and an intermediate expansion is required. This region is of width 
0(e), and the match can be achieved if we choose 

e - 1/ inU"*) . (24) 

The variables to be used in this region are defined by 

y o cin X , Q(y) * H/X , (25) 

and the governing equation (20) becomes, after one integration, 

a + Q(Q + e'y/') (c3 iifi - E |2) = 0 . (26) 

3y 3 <jy 

Neglecting the exponentially small term, the solution Is 

» Cq + cc^ + 3(Uj^ + CU 2 ) y + O(c^) , (27) 

which matches with the inner solution (23) if we choose c^ » 1 , Cj = 0 . To match with 
the outer solution (18) we write 

y=l+cln(a-x) (28) 


and obtain 


h 'X' (a - x) (1 + + e(1 + 3u^)“^^^(Uj ln(a - x) + U 2 ) + O(c^) 


Comparing this with the outer solution, we can see that the logarithmic term and the 
constant term both match if 


1 + 3u, 


, U 2 = u^ (2 - In 2a) , 


so that the rate of spread of the drop is given approximately by 
a ■ c{l + c(2 - in 2a)) (a*^ • 1) , 
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\ 


i 


or, to the same order of accuracy, 


3 



da 

at 


a 


-6 


1 . 


(32) 


The objective of determining the rate of spread of the drop as a function of its width 
has thus been achieved. If the leading term only in the outer expansion is included an 
apparent contact angle is given by the value of a"2^ from (18). With variables in their 
dimensional forms, equation (32) yields the expression 


a 


a 


{■ 




(33) 


for the apparent contact angle at the edge of a drop of width 2a advancing with speed U 
when the static contact angle is small. But the inclusion of the second term in the outer 
expansion destroys the validity of the concept of an apparent contact angle as a directly 
measurable quantity since, with the extra term included, the slope does not tend to a 
constant at the edge of the drop within the outer region. It is, however, still possible 
to define the quantity (33) as a derived contact angle as it could be measured by finding 
the curvature of the drop at its mid-point, for example, thus taking into account the outer 
solution but not requiring any measurement to be made in the vicinity of the edge. 


Sliding 

The second problem is when the drop is placed on an inclined plane. The drop is no 
longer symmetric and both edges have to be treated. As before, the Bond number is small, 
and gravity enters the problem only through the positive parameter K defined in (4) , 
which measures the component of gravity down the plane. Any transient behaviour is again 
ignored, so that we suppose the height of the drop to depend on the co-ordinate x 
measured along the plane and on the positions of the two edges, but not directly on the 
time. The same three regions encountered in the spreading problem are present, but 
because both edges have to be treated, it is convenient to introduce a change of variable. 
If we write h(x,t) = f(s,a^,a 2 ), where 

X = ^(a^ + a 2 ) + ^(a^ - a 2 ) s , (34) 

the range covered by the independent variable s is from -1 to 1. In the outer region, 
we c \T\ write f = f^ + f^, + ... , and the equation for f^ is 


T + K = 0 , (35) 


and the solution which vanishes at the edges and which satisfies the volume condition (6) is 


f 


0 






(36) 


Since f^ must be non-negative for s in [-l»l] , this solution is only acceptable if 
0 < b i 1 , where 

” TO ■ *2’"* 


This condition may be broken either in the initial stage or during the spreading of the 
drop. In either case, a different appr*ach from that used here must then be employed. 

The failure of the condition implies that regions are developing where the drop becomes 
very thin and although this is an interesting possible behaviour, this aspect of the prob- 
lem is not investigated here, and we assume that b < 1 throughout the motion. 


The second term in the expansion in the outer region satisfies the equation 

a^d+sMltb) + 42(1-8) (1-b) 

(aj-aj)* a-s2)2(i + bs)* 

The be .ndary conditions on f ^ are 

fj * 0 at s « -1 and at s«*l, |fjds*0. 


(38) 


(39) 
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and the solution can be found In closed form after a great deal of algebra. The only 
quantities of interest are the asymptotic values for the height of the drop near the two 
edges and these are given by 


4(a^-x) (1+b) 


(a^-a 2 ) ** (a^-x) 


/ (l-b)(a.-x) 


; ‘ 1 


(a,-a2)^ 

X 

16(l+b)^ 


(1+b) (a^-aj) 

+ 

T=F| 


t 

4(x-a2) (1-b) 



- 

( (l+b)(x-a,) 

, l+2b 1 

+ i -i- 1 


(a^-aj) 


16(l-b) ^ 

"®2 

(1-b) (a^-a2) 


®i T+E 



(40) 


(41) 


The solutions in the inner and intermediate regions are similar to those found in the 
spreading problem, except that the edges may be moving in either direction or be station* 
ary. The conditions (7) apply at the lower edge x * aj and when they are applied, and 
the resulting inner expansions matched via the intermediate expansion to the outer solu- 
tions, we obtain the equation 



(1+b) 

xTl-b) 


S(a^) 


i-2b 

-T=F 


\ 

1 • ®2 

1 r2®(i+b)^ - (A ) 1 

/ 

I 1 , ,6 ® '®1^ f 


(42) 


where S is a step function, defined by 

S(x) = 1 for X > 0, S(x) = B for x < 0 , (43) 

When the edge is at rest, (42) is replaced by » 0 . The corresponding results for the 
upper edge x = , where the conditions (8) apply, are 


in 



-a^) (1-b) 
X(l+b) 


S(-a2) 


l+2b 

1+ET 



i+b 



2Ml-b) 

(aj-d2) 


3 

T 




(44) 


when the edge is moving, and ^ 0 when it is at rest. 

These equations are sufficient to determine the future behaviour of the drop from any 
given initial position of the edges. The drop may not move at all, or one or both edges 
may move in either direction as the drop spreads or contracts. The most interesting 
possibility is when the drop spreads and slides, approaching a final state in which the 
width of the drop and its ^ speed down the plane attain constant values. This final state 
can be found by setting a^ = *= U and a.-aj « 2a . From (42) and (44) , with b re- 

placed by its value in terms of K and a irom (37) , we obtain the equations 


U * 


■^(I + a^K'^) - 1 + 

r.n/1^) . 

\ \ / 


12a°K^ 

9 - a®K^ 


(45) 


* 3-a-K 


, /2a 3-a^K 

f— - — r- 
X 3+a'*K 


J 


6a‘*K 

3-a*' -^ 

6a**K 

3+a**K 


a^ - (1 + ya**K) 


(46) 


(1 


3 a**K) 




There are too many parameters for it to be easy to make general statements about the 
motion. An approximate set of criteria can be found by retaining only the dominant terms 
on the left-hand sides of (42) and (44) , that is, the logarithmic terms. Then it follows 
that 

a. 

(47) 


‘1 


if 

if 


B < a~2 (1 + a'*K/3) < I , 
•■*2(1 - a**K/3) < 1 . 


< a 
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A static final state is only possible if these two conditions on the width of the drop 
overlap, which they do if 0 « X ^ ^ where * 

*c " 

while for K > the drop will slide. If there is no contact angle hysteresis, that is, 
if Of * , the drop will slide however small the inclination of the plane • 

Extensions and conclusions 

Although only thin two-dimensional drops have been considered here, the method is 
applicable to many other cases. The spreading of a thin drop by both capillarity and 
gravity and without the restriction to two-dimensionality has been examined**. The 
capillary spreading of a drop which is not thin, so that the simplifications of lubrication 
theory are not available, has also been examined and the results compared with those 
obtained experimentally^. Further work on the sliding problem, with the extension to 
three-dimensio Aality and the lifting of the restriction on the gravity paraiaeter K is 
planned. 

The aim of this paper has been to show that the moving contact line problem can, in 
certain circumstances, be solved in a satisfactory manner. The solutions for which 
experimental corroboration is available^ indicate that the proposed boundary conditions 
can be used with some confidence. The results obtained here and elsewhere** ^ show that the 
application of these conditions need not involve a large amount of refined numerical 
analysis to resolve the solution near the contact lijie. 
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